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Abstract 

The slope of the Isgur-Wise function at the normalization point, £^(1), is one of the basic param- 
eters for the extraction of the CKM matrix element V c b from exclusive semileptonic decay data. 
A method for measuring this parameter on the lattice is the effective theory for heavy quarks at 
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small velocity v. This theory is a variant of the heavy quark effective theory in which the motion 
of the quark is treated as a perturbation. In this work we study the lattice renormalization of 
the slow heavy quark effective theory. We show that the renormalization of ^^(l) is not affected 
by ultraviolet power divergences, implying no need of difficult non-perturbative subtractions. A 



lattice computation of ^'(l) with this method is therefore feasible in principle. The one-loop 
renormalization constants of the effective theory for slow heavy quarks are computed to order v 2 

•rH , 

^ ■ together with the lattice-continuum renormalization constant of £^(1) • 

We demonstrate that the expansion in the heavy-quark velocity reproduces correctly the in- 
frared structure of the original (non-expanded) theory to every order. We compute also the one- 
loop renormalization constants of the slow heavy quark effective theory to higher orders in v 2 
and the lattice-continuum renormalization constants of the higher derivatives of the £ function. 
Unfortunately, the renormalization constants of the higher derivatives are affected by ultraviolet 
power divergences, implying the necessity of numerical non-perturbative subtractions. The lattice 
computation of higher derivatives of the Isgur-Wise function seems therefore problematic. 
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1 Introduction 



The effective theory for heavy quarks (HQET) []l], (for a comprehensive review and refer- 
ences to the original literature see ref.||) allows a clean determination of the Cabibbo-Kobayashi- 
Maskawa matrix element | V c b | from the exclusive semileptonic decays of B mesons 

B -> L> w +l + ui (1) 

recently measured by the ARGUS ||] and CLEO || collaborations. 
In the limit of infinite mass for the charm and the beauty quark, 

m c , m b -> oo, (2) 

the six form factors parametrizing the hadronic matrix elements of the decays (Q) can all be ex- 
pressed in terms of a unique form factor, the Isgur-Wise function £(v ■ v') || 0, ||, 

(D,v\jfr c (0)\B,v'} = VM D M B ( Vfi + ^) i{v ■ v') (3) 
(D*,v, e | J^ c (0) | B, v') = -y/M D M B [ie^^v'^ + e M (l +vv')- v^v' ■ e)({v • v') (4) 

where v' and v denote respectively the b and c quark 4-velocities and jj^ c (x) is the weak current 
describing the transition of a beauty quark into a charm quark, J b ^ c (x) = c(x)7 M (l — 75)6(3;). 

This function is normalized at zero recoil 

t(v'-v = l) = 1, (5) 

and 1/m-corrections vanish in this kinematical point |J. A model independent analysis of the 
decays (|l]) extrapolates the experimental data up to the endpoint, where the form factors are 
known by symmetry. In order to eliminate the systematic errors introduced by the extrapolation, 
it is essential to know also the derivative of the Isgur-Wise function at the normalization point, 
£«(!). 

So far, four methods have been proposed to compute on the lattice the slope of the Isgur-Wise 
function. The first one was the estimation of the derivatives of the Isgur-Wise function at the zero- 
recoil point from the lattice determination of the Isgur-Wise function at discrete points in a region 
close to the zero-recoil point flTCfl . This method presents some problems related to the extrapolation 
to the zero-recoil point that could lead to uncertainties in the determination of ^^(l). The authors 
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of |Tl| have suggested a new method to compute directly on the lattice the slope of the Isgur-Wise 
function which avoids any kind of extrapolation. They proposed to study the first spatial moments 
of two- and three-point meson correlators. Then the derivatives of the Isgur-Wise function could 
be extracted by forming appropriate ratios of these correlators. The UKQCD Collaboration has 
recently carried out a exploratory study of the feasibility of this method [^] . The main conclusion is 
that there are large finite- volume effects in the lattice evaluation of the moments of the correlation 
functions, having a geometrical origin. Therefore, by increasing the length of the lattice in the 
spatial directions, these undesiderable volume effects can be reduced. Unfortunately, they are large 
on currently available lattices. Some approximations have been presented by this group in order to 
control the volume effects and extract the slope on finite volumes. 

Both computations described above treat the heavy quark as an ordinary quark but with a 
small hopping constant. The first calculation of the Isgur-Wise function using the lattice HQET 
has been done by the authors of O] . In that work, the lattice propagator of the heavy quark with 



velocity v is obtained from a Wick rotated lagrangian [14]. 

The fourth method is based on an expansion of the HQET in the heavy-quark velocity around 
the static theory (hereafter called 'effective theory for slow heavy quarks', SHQET) 17]. As in the 
method of spatial moments the derivatives of the Isgur-Wise function at the zero-recoil point 
can be extracted directly from ratios of two- and three-point correlation functions (see Section 2 for 
details). The main point is that there are no unexpected geometrical volume effects in the lattice 
computation of these correlators because the static propagator is local in space. Moreover, the 
SHQET circumvents the problem of the euclidean continuation of the Georgi theory for heavy- 



quarks with non-vanishing velocity [14, 15|. In fact, in order to simulate heavy-quark with velocity v 
on the lattice, the continuum Minkowskian HQET must be transformed into a discretized euclidean 
field theory. The analytical continuation is not simple because the energy spectrum is unbounded 
from below [14, [16|]. On the contrary, expanding around small velocities, we are perturbing the 
static theory whose energy spectrum is bounded from below. Roughly speaking, we may say that 
the heavy quark has a 'perturbative motion' in the SHQET produced by the 'velocity operator' 
(v ■ D). This theory has not been used yet in numerical lattice simulations. 

In this paper, we analyse the lattice version of the SHQET. One of our main results is that the 
lattice renormalization constant of £^(1) does not contain any ultraviolet power divergence (i.e. 
proportional to l/a n , where a is the lattice spacing). The renormalization of £^(1) is affected only 
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by logarithmic divergences (of the form log am) which can be subtracted with ordinary perturbative 
computations. This implies that the lattice computation of £^ (1) by simulations using the SHQET 
is feasible in principle. Ultraviolet power divergences are indeed a serious problem for numerical 



simulations, because they cannot be subtracted perturbatively with adequate accuracy |18|, [19]. 
We also compute the one- loop lattice renormalization constant of £W(1). The knowledge of this 
renormalization constant is essential for converting the values of £^(1) computed on the lattice to 
the values in the original (high-energy) theory. Moreover it is shown that the infrared as well as the 
ultraviolet behaviour of the non-expanded theory (HQET) are reproduced by the SHQET order 
by order in the velocity. This is a non trivial check of the consistency of our approach. We also 
give the lattice renormalization constants of higher derivatives of the Isgur-Wise function, £( n )(l), 
n > 1. Unfortunately, the lattice renormalization constants of higher derivatives £( n )(l) for n > 1 
are affected by power divergences. The computation of higher derivatives with SHQET is therefore 
more difficult than that one of ^^(l). 

This paper is organized as follows. In section 2 we review the SHQET in the continuum 



[17]; the derivatives of the Isgur-Wise function are expressed as ratios of three- and two-point 
correlation functions. Section 3 deals with the lattice regularization of the SHQET. In section 

4 we briefly review the matching theory of lattice operators onto the continuum ones. In section 

5 we renormalize the lattice SHQET at order a s and to all orders in the velocity. In section 6 
the renormalization of the heavy quark current Jj^° is computed. In section 7 we calculate the 
lattice-continuum renormalization constants for the derivatives of the Isgur-Wise function. Section 
8 deals with the problem of power divergences. Finally, in section 9 we present our conclusions. 
There are also two appendices where the technique to compute lattice integrals is described. In 
appendix A the analytical expressions and numerical values of one-loop diagrams are presented and 
in appendix B the method for subtracting infrared divergences is briefly explained. 

2 The effective theory for slow heavy quarks 

In this section we review the basic results and formulas of the SHQET. The Georgi lagrangian 



describing a heavy quark Q with velocity = (vT-j- v 2 , v) in Minkowsky space [0] 

C(x) = Q\x)iv ■ D(x)Q(x) (6) 



is decomposed as 



£(x) = Cq(x) + Ci(x) 
3 



(7) 



where 

C (x) = Q\x)iD (x)Q(x) (8) 
is the static unperturbed lagrangian and 

d{x) = Q\x)i[D (v - 1) - v • D}Q(x) (9) 
is a perturbation lagrangian giving rise to the motion of Q. 

From this splitting it is easy to derive the following expansion of the propagator of Q [|l7]] 

S(x,y;v) = -iG(t x -t y ) ( P(t x ,t y ) + /*" dt z P{t x ,t z )v ■ D(x,t z )P(t z ,t y ) 

\ Jty 

rtx rt z _ 
+ / dt z dt w P(t x ,t z )v ■ D(x,t z )P(t z ,t w )v ■ D(x,t w )P(t w ,t y ) 

J ty tJ ty 

-f-P(t x ,ty)+...jS(x-if) (10) 

where P(i&, t a ) is a P-line in the time direction joining the point (x, t a ) with the point (x, tb) 

P(t b ,t a ) = Pexp^ig A (x,s)dsJ (11) 

S(x, y; v) is expressed as a sum of static propagators with an increasing number of local insertions 
of (v ■ D) giving rise to the 'perturbative' motion of Q. 

The spin structure of Q is taken into account multiplying S(x,y;v) by (1 + v)/2 



1+ / 

H(v) = —t- S(v) 



1 + 7Q o( ) , ( 1 + 70 5 (i) _ J3 si0) \ 
2 \ 2 2 J 



"3 



+ (^SV _ + | 5 (0) j ,2 + (,1) (12) 

where we have taken the heavy-quark moving along the z-axis. Inserting the propagator H(v) in 
the Green's functions describing the dynamics of heavy flavored hadrons, we have the following 
expansion in powers of v% 

G(v) = G^ + G^ v 3 + G^vl + --- (13) 
Consider now the following three- and two-point correlation functions 

C 3 (M') = jd 3 xd 3 x'(0\ T [oUx'),jfr c (x),O B (0)] | 0) (14) 
C B (t) = J d 3 x(0 | T [O B (x),O B (0)} | 0) (15) 
C D (t' - t) = j d 3 x'(0 | T [tf D {x'),0 D {x)\ | 0) (16) 



where Oh(x) is an interpolating field for the H meson. The simplest choice (which we adopt in 
the following) is: Oh(x) = Q(x)i'y5q(x), where Q(x) = b(x), c{x) for B and D mesons respectively, 
and q(x) is a light quark field. For large euclidean times, t — ► oo, t' — t — * oo, the matrix element 
is given by 

C 3 (M') 



{D,v\J<r°{0)\B,v') 



C B (t) C D {t'-t) 

where Zb and Zb are the renormalization constants of the operators Ob(x) and Ob(x), given by 



^ = (B,v' | O B (0) I 0) 
%^ = (D,v | O D (0) | 0) 



(18) 



Since both the wave functions and the interpolating fields in eqs.(^) are pseudoscalars, the matrix 
elements do not depend on the velocity v(v'), unless we deal with smeared currents [jD]. 

Inserting the propagator (|l2[) for the c quark in C^(t,t') and Co(f — t) we derive expansions 
of the form (|l3|). Inserting them into eq.(17), we get the following expression for the derivatives of 
the Isgur-Wise function with respect to V4 at the zero recoil point V4 = 1 
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where we have used the identity 
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j« <2) (i) 
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(19) 
(20) 



(21) 



Higher derivatives can be computed similarly. 



3 Lattice regularization 

We consider the discretization of the effective theory for heavy quarks proposed in ref.|14], for- 



ward in time and symmetric in space. For a motion of Q with velocity along the z-axis v 



(0, 0, Vs, J 1 + uf), the action iS is given by 

iS = - v^{x) [ip(x) - Ul(x)ip{x - 4)] + 

x 



-i^(x) [U 3 (x + 3)^(x + 3) - C/ 3 T (x)^(x - 3)] (22) 



where fx is a unit vector in the direction li, and U^{x) are the links related to the gauge field by 
Un(x) = exp[-igAfj,(x - fx/2)]. 

The Feynman rules are those of the static theory plus additional interactions generating the 
motion of Q. Assuming a convention for the Fourier transform according to which ip{x) ~ ex.p(ik-x), 
we have 

^ (O) (*0 = i hkT— (23) 

1 — e lKi + e 

^ 0) = ig 6» 4 t a e-^+ k '^ (24) 
V$ )tod = -^^4*0*6 e" ifc * (25) 

The linear interactions in v 3 are given by 

V (1) = -v 3 sin fe 3 (26) 
V™ = gv 3 6^t a cos(k 3 + k' 3 )/2 (27) 
Vil )tad = ^S^ 3 S u3 t a t b smk 3 (28) 

and the linear interactions in («4 — 1) are given by 

y (2) = -(v 4 -i)(i- e - ifc4 ) (29) 

^ = ^(^4-l)^4^e-^ 4+fc ' 4 )/ 2 (30) 
V$ tad = -^-(v 4 -l)S^S v4 t a t b e- ik ^ (31) 

where k and fc' denote respectively the momenta of the incoming and outgoing heavy quark, and 
is the interaction vertex of the heavy quark with a gluon provided with a polarization along 
the fx axis. Vj^f" are the vertices for the emission of two gluons, for the case of the tadpole graph 
(k = k'). We notice that the vertices labelled contain second and higher orders in the velocity 
v 3 . It is convenient to keep them unexpanded. Finally, note that the conventions for the sign of the 
Fourier transform and of the velocity are not independent, if one wishes to intend k as the residual 
momentum of the heavy quark. 



4 Renormalization of lattice operators 

Since the lattice effective theory and the continuum one are two different versions of the same 
physical theory, the matrix elements computed in both theories must coincide. This is a non-trivial 
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condition to impose (matching condition). We match amplitudes of the bare lattice theory onto 
the corresponding ones of the continuum theory renormalized in some chosen scheme (such as for 



example MS). If the lattice lagrangian and the continuum one have at the beginning the same 
parameters (masses, couplings, etc.), matching is accomplished adding appropriate counterterms 
to the lattice lagrangian. If we are interested also in the matrix elements of composite operators, 
an analogous matching has to be performed: appropriate counterterms have to be added to the 
lattice composite operators. Because of mixing, to a renormalized operator in the continuum it 
corresponds in general a linear combination of lattice bare operators. 

The technique for obtaining the lattice counterpart of a continuum operator is standard |p0| , 
|2l"|| . In lattice regularization the inverse of the lattice spacing 1/a acts as an ultraviolet cut-off, 
and bare lattice amplitudes depend explicitely on a. Continuum amplitudes depend instead on a 
renormalization point /u. To avoid large logarithms in the matching constants (loga/x >> 1), let 
us first match the amplitudes by taking fj, = a . At this stage we are therefore dealing with the 
finite discrepancies coming from the use of different regulators. The relation between continuum 
and lattice operators at one-loop level for fi = a -1 is given by 



a.Ja ] 



= a" 1 ) = £[<% + 5Z l3 ] Of *(a) (32) 



7T 

3 \ 



where Of ont are the operators in the continuum we are interested in, 0^ att are the lattice ones 
and 5Zij are finite renormalization (or mixing) constants. The sum extend over all the operators 
that can mix with 0[ att as a consequence of the symmetry breaking induced by the continuum 
and the lattice regularization procedure. By sandwiching the operators between arbitrary external 
states of momenta p, we derive the following relation involving Green's functions of the bare lattice 
operators or the (renormalized) continuum ones 

Usia- 1 ^ ^Cont,Latt {p)] < Q . >(0 ) (33) 



< o c r tMtt >= £[<% + (^p) eg 



J 

where the superscript (0) denotes tree level matrix elements. Demanding compatibility between 
d32| ) and d33|), we derive 

SZa = ]im [c^ n \p) - C^ u (p)} (34) 

The mixing coefficients 5Zij are independent of both the external states and the momentum con- 
figuration used to calculate matrix elements of 0%. 

We consider now the matching in the more general case /ia/ 1. Since we already matched the 
amplitudes at \xa = 1, we need only to evolve the mixing coefficients 5Zij from /i = 1/a to a generic 
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renormalization point with RG techniques. At one loop-level, the SZij's do depend on the renor- 
malization scheme. In order to obtain a renormalization scheme independent matching condition, 
the two-loop anomalous dimension contribution must be taken into account in the diagonal terms 



22 



n Cont (lA ( »s(a r ) Y l/(31 ( a s {a *) a s (v) \ , 

0< {fl) = \^w) [ + {~ —) RoA 

x + (^P) ^\O h 3 at \a) (35) 



j 



where 



ROi = 4 [72 A -71/%] (36) 
Pi 

with j n the n-loop anomalous dimension of the operator Oj defined by 
Of = Z Q Of 

7 = -V-^logZo = 71 (^) + 72 (^) + ••• (37) 
and with (3\ and the one-loop and two-loop coefficients of the /3-function respectively 

(3(a s ) = + + ... (38) 



We have 



11 1 

A = "T + 3 nF 

51 19 

ft = -T + i2 nF (39) 



The expression for the running coupling constant is given by 

2(3 2 log log^/A 2 ) 



2^ 



A log(/^/A 2 



1 + 



ft 2 log( M 2 /A^ 



(40) 



where np is the number of active quark flavors and we can take A = 200 MeV in the MS scheme. 

It is expected that for the values of a -1 currently used in lattice simulations the matching 
should depend only weakly on the continuum regularization. This is why physicists usually compute 
matching constants without including the two-loop anomalous dimension. 

Finally, let us briefly expose the problem of the power divergences in lattice computations. Since 
QCD is asymptotically free, the matching constants SZij can in principle be computed with RG- 
improved perturbation theory in the limit a — > (in practise one requires aA <C 1). Unfortunately, 



the mixing coefficients contain in some cases inverse powers of a, which divergence as a goes to zero. 
Then, in computations of matrix elements of the continuum operator, the leading term is this mixing 
term of 0(l/a n ) which is a lattice artifact generated by the regularization procedure and thus must 
be subtracted. In other words, in order to obtain finite Green functions of composite operators on 
the lattice, we must subtract power divergences in a -1 from the Monte Carlo data. It has been 



argued that it can be done in perturbation theory. However, as pointed out in refs.fllq, 19, ^4|, 25], 
it is not clear that the coefficients of power divergences can be calculated to sufficient accuracy 
in perturbation theory. In general, very difficult non-perturbative subtractions for lattice Green's 
functions are required. 

5 Renormalization of the lattice SHQET 

In this section we discuss the renormalization of the SHQET given by the lagrangian (|22"|), i.e. the 
determination of the counterterms which have to be introduced to match amplitudes of the lattice 
SHQET onto those of the continuum HQET. 

To obtain the renormalized operator (v ■ D) we compute the one-loop heavy quark self-energy 
with insertions of (v-D) using the lattice Feynman rules of section 3. This is equivalent to calculate 
the one- loop heavy quark self-energy up to a given order in the velocity v^. 

After that, we match the heavy quark propagator of the lattice SHQET onto the continuum 
HQET propagator, expanded in V3. 

For calculational convenience, we will take equal incoming and outgoing momenta and adopt the 
Feynman gauge for the gluon propagator. The infrared divergences which appear at zero external 
momenta are regulated giving the gluon a fictitious mass A. No problem arises with non-abelian 
gauge symmetry because all the amplitudes are QED-like. Other choices are possible for the 
infrared regulator, such as for example to take virtual external states [23]. However, by using a 



non-vanishing gluon mass, we achieve a great simplification in computing the lattice loop integrals. 
Indeed, we can safely Taylor expand the corresponding diagrams about zero external momenta up 
to order 0{a) to determine all nonvanishing terms as a goes to zero. Upon doing this, we will 
subtract the infrared (logarithmic) divergences from the integrals with the technique explained in 
Appendix B. 

The computation of the diagrams will be done with two different methods for dealing with 
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the non-covariant poles coming from static lines. The first method is based on partial integration 
with respect to k<± (the fourth-component of the euclidean loop momentum) in order to eliminate 
the non-covariant poles. The integrand is reduced to a covariant form and can be computed with 
usual techniques. The second method is to integrate analitically over k& using the e-prescription of 
the static heavy-quark propagator and the Cauchy's theorem. This latter technique involves less 
algebra, but leads to non-covariant 3-dimensional integrals. The comparison of the results obtained 
with the two methods provides us with a check not only of our analytical computation but also of 
our numerical calculations. 

5.1 Heavy Quark Self-Energy up to 0(i>§) 

To illustrate the method of partial integration, in this section we briefly describe the computation 
of the diagrams that determine the heavy-quark self-energy up to 0(v 2 ), which are depicted in 
Fig.l. We do not consider the insertion of ^-vertices (i.e. those in eqs.(p9|) to (|3l|)) because their 
contribution can be shown to be trivial. We will treat this subject in detail in the next section. 

We start by computing the diagrams with one insertion of the operator (v-D) in Fig.l. Diagrams 
A.l and A. 2 vanish in the Feynman gauge. This happens because the gluon is emitted by the 
operator {v ■ D) with a polarization along the z axis, while it is absorbed by the static vertex with 
a polarization along the time axis. The (amputated) amplitude of diagram A.3 is given by 



where p is the external momentum, Cf = ^2t a t a = (N 2 — 1)/2N for an SU(N) gauge theory, and 



because the integrand is odd in k 3 . 

First derivatives of A${p) with respect to the external momentum contain logarithmic ultraviolet 
divergences. The only non- vanishing derivative is that one with respect to p%. With a partial 
integration with respect to k^ of the factor 




(41) 



Ai(0 = E^=il-cos/ M + (aA) 2 /2. 



Now, Az(p) vanishes at zero external momentum p = 0, 



A 3 {p = 0) = 0, 



(42) 



— ik4 



(1 - e" ifc 4 + e )2 



(43) 
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we reduce the integral to the following form 



dps/ o 167r 2 6tt 2 J-k Ai(k) ; 
where r/(A;) = X)f=i srn2 ^« an d a symmetrization over the spatial momenta has been done. The 



infrared singularity of the integral is isolated with the technique introduced in ref.[23] and described 
in detail in ref.|l6|]. The result can be written as 

'8A 3 ' 



^) ^f V3 [-2log(aX) 2 + a 3 ] (45) 

7T / 4 



9P3 /o 

where the subleading (finite) term 03 is a constant evaluated numerically, a 3 = 0.448. 
The tadpole graph A. 4 is given by 

, f s 9 2 C F vs . d 4 k 1 

Mp) = ~^ smp3 L (^fatm (46) 

As in the case of diagram A3, A^(p) vanishes at zero external momentum, 

A 4 (p = 0) = 0. (47) 

The first derivative with respect to p 3 is finite (i.e. does not contain logarithmic divergences) and 
reads 

dA 4 \ ( a s \ C F 

— — v 3 a 4 (48) 



dp 3 ) W / 4 
where 04 is a numerical constant, 04 = 12.23. 

Let us consider now the renormalization of a double insertion of (v ■ D) at zero momentum. We 
study the Green function 

G(z, w) = { d 4 xd A y (0 I T[Q(z) Q\x)(v ■ D)(x)Q(x) Q\y)(v ■ D){y)Q{y) Q\w)} | 0) (49) 



The diagrams involved are drawn in Fig. 2. Diagrams B.l and B.2 vanish in the Feynman gauge. 
The amplitude of diagram B.3 can be written as 

9 2 C F vj _ ipi /■+* d 4 k fj(jfe) e- ik * 



B3(p, = -°—l e -n L (50) 

The amplitude at zero external momentum is given by 

fl 8 (0) = ,|4 r ^ 7 f } C X n ^ ( 51 ) 

iy ' 16^ 2 3 6TT 2 J- n (l-e- lk 4 + e) 3 Ai(k) v ; 

It is convenient to reduce the integrand to a covariant form, by eliminating the triple pole coming 

from the static line. Let us describe in detail the transformation of this integral, which will illustrate 

the technique to deal with poles of odd order. 
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We perform first a partial integration with respect to k^ analogous to that one of dA^/dp^. 
This transformation brings the integral into the form 

B(ti\ 9 2 C F 2 -l [+* 4 r?(k) (e** + 1) 1 

B3{0) ~ W^ V *2^L k (l-e-* 4 + e) Al (k)2 (52) 
The simple non-covariant pole is treated by writing |2 

1 ( 1 1 \ 1 

+ ^T^T (53) 



Aiik) 2 yAi(k) 2 Ai(0,fe)V A!(0,fc) 2 
In the integral containing the difference of gluon propagators, 

/+7T 1 I p ik4 / 1 1 \ 
d 4 k r?(k) — — ttttt , (54) 

l_ e -ik 4 + e 'V 7 ^ Al(k )2 Al(0)k)2y /' v > 

one can set e = 0. Since the gluon propagator is even with respect to k%, one can symmetrize the 
factor 

■ 1 + cos &4 (55) 



The integral / therefore reads 



i = r d^ (i+c A ° s ^ )??(g) - 2. r d 3 k (56) 

i-Tr Ai(k) J-K Ai(0,k) 

In the remaining integral 



J = r d k 4 1 + f Z+Vk ^ , (57) 
i-Tr 4 l-e" lk 4+e Ai(0,k) 2 

we perform the contour integration over fc 4 analytically by setting z = exp(ik^). 1?3(0) is finally 
expressed as a sum of a 4-dimensional integral and a 3-dimensional one 

16vr 2 3 ^ 24vr 2 7 A^) 2 12W A^O,*;) 2 / 



The integrals in eq.(|58|) are infrared finite and are easily computed numerically 

Bs(0) = (f) ^ «I &30 (59) 

where 630 = —5.044. 

The first derivative of B^(p) with respect to p^ is logarithmically divergent, and is given by 

dB 3 \ „ . , g 2 C F v\ f d A k smk 4 e~ ik * rjCk) , , 



Op 4 /o " 6 J (2f) 4 (l-r ,k H£) 3 Ai(fc) 2 
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Using the same tricks as for .63(0), this integral is transformed into 



dEh 
dp a Jo 



-iB 3 (0) + 



+ 



9 2 Cf 2 i 



127T 2 



/ d 4 k [1 + 2 cos k 4 + cos 2k 4 

V(k)_\ 
kf J 



(61) 



12vr 2 7 A!(fc) 2L 7 4J 12W Ai(0,< 

The logarithmic divergence of the amplitude (the log(aA) term) is entirely contained in the first 
integral. The computation yields 



(dB 3 \ . fa s \ C F 2 



21og(aA) 2 + 6 3 i 



(62) 



where 631 = 4.988. 



Finally, the amplitude of diagram B4 is 

g 2 C F vl [+* d 4 k 3 + a(k) 



B 4 (p) = 



12 



/• +7r d"k 

7_ w (^o 1 r 



e -ifc 4 + e Ai(fc-p) 
where a(k) = J2i=i cos h- 

The computation of -B 4 (p) is analogous to that of B 3 (p). We have 



(63) 



y ' 16vr 2 ^247r 2 J A x (k) 12vr J Ai(0,k)/ 



Upon a numerical computation we find 

54(0) = (^) 6* 

where 640 = 20.566. 



(65) 



The derivative with respect to p 4 reads 



dB A 



9 2 C F 



-i f +7T , 

24tt 2 y_ w 



dpWo !6vr 2 3 24vr 2 
and the corresponding numerical computation yields 



+ "^4 k [3 + a(k) ][ l + cosk 4 ] 



Ai(k) 2 



(66) 



(67) 



where &41 = -2.485. 

Putting all contributions together, we can write the heavy-quark self-energy up to order 0(v 3 ) 



as 



S(p,u) 



1 /q s \ Cjr 



a V 7T 



y(°) + „ 2 y( 2 ) 

ZJQ -f- f 3 Z^Q 



«4 



+ 





C F r 


(?) 


4 




C F r 


(?) 


4 - 



Sfj + « 3 2 sg - 4 log(aA)] (ip 4 v A ) 
E$ - 4 1og(oA)l (p3«3) 



+ 0(a s a, fg 



(68) 
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where £ °^ and S^q are the mass and wave function renormalization of a static heavy quark [26] 



Their numerical values are tabulated in Table A.l. On the other hand, the numerical values of the 



new constants Eg^o 40 are 



4 2) = ho + ^40 = 15.52 
S40 = hi ~ ho + hi = 7.55 



£<J = a 3 + 04 = 12.68 (69) 



In the next section we compare (69) with the self-energy calculated by using a different integration 
technique. 

5.2 Heavy Quark Self-Energy beyond 0(v%) 

Here we compute the heavy-quark self-energy at one-loop in the coupling constant a s but to all 



orders in the velocity V3. We will demonstrate that the discretized lagrangian (22) reproduces the 
correct infrared behaviour of the HQET to all orders in the velocity, as it should be. The calculation 
will be performed utilizing a different method from the one used in the previous section. This is 
useful to check both our analytical and numerical results. 

We start by noting that at one-loop the diagrams that contribute to the self-energy of the 
heavy quark at order 0(v™) are those depicted in Fig. 3. In fact, these diagrams represent the only 
two ways of inserting m t^-vertices (i.e. those in eqs.(p6|) to (p8|)) into the gluon-loop self-energy 
diagram. The reader may however argue that we are ignoring the considerable number of 1)4- vertex 
insertions (i.e. those in eqs.(p9|) to (^TJ)) which give rise to corrections to the heavy-quark self-energy 
of the same order in the velocity as those considered above (see Fig. 4). This is of course true, but 
it is very easy to show that the full effect of all possible ^-vertex insertions into the diagrams of 
Fig. 3 is just to multiply them by a factor (l/tti) 771 . 

To demonstrate this result, consider the diagrams of Fig. 3. We can insert n ^-vertices only 
in two ways, namely, A: one at any of the quark-gluon vertices and the remaining n — 1 on the 
heavy-quark line and B: all n insertions into the heavy quark propagator inside the gluon loop. 
The resulting diagrams are shown in Fig. 4. Now, the effect of a ^-vertex insertion into a heavy 
quark propagator is simply to multiply the same propagator by (— ) (1)4 — 1), as it can be seen from 
the Feynman rule in (^). If the insertion is at a quark-gluon vertex, the effect is to multiply the 



same vertex by (v^ — 1) (see eq.(|30D). For example, two ^-vertex insertions, one at the quark-gluon 
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vertex and the other into the heavy-quark propagator inside the gluon loop, give (— ) (v± — l) 2 times 
the old diagram without any t>4-vertex insertions. Therefore, we only have to count the number of 
topologically different diagrams with n t>4-vertex insertions in each class of Fig.4, for all diagrams 
in a class give the same contribution to the heavy-quark self-energy. To this end, we observe that 
the number of different ways we can insert n t>4- vertices on a heavy-quark line where there are m 
i>3-vertices is (n + m)\/[n\ ml] (the old combinatorial problem of distributing n balls inm+1 boxes). 
Therefore, the sum of the graphs in Fig.4 gives 

Fig.4 = (-)" {" ± m " 2 | ! (v 4 - IT C.l + (-)" (V4 ~ IT C2 

[n\ (m — 2)!J [n\ m\\ 

+ 2 { ^~ l t m) l («4 - IT C2 + 2 { H~ 2 t m) l (v, - IT C2 

v ' [(n-l)!m!] v 1 v ; [(n-2)!m!] v ; 

= { - )n mZ-m (Vi - lT(CA + c - 2) m 

Summing from n = to oo, we get 

]T Fig.4 = — (C.1 + C.2) (71) 

n=0 V 4 

as anticipated. 

We turn now to the computation of diagrams of Fig.3. Their amplitudes are 



dz z —z 



2m z [1 + e - e~ tp4 z] m+l (z - z-)(z - z+) 



(72) 



and 



C 2 (m,p) = - C F (-v 3 r±- r sin m - 2 (p3-fc3)cos 2 (p 3 -A;3/2) 

^ (73) 



2vri z [1 + e - e' ip4 z\ m - x (z - z-){z - z+) 

where z = e lki , the contour integral is along the unit circle and we have used the fact that the 
gluon propagator can be written as 

1 -z 



Ej=i (1 " cos(^)) + (aA) 2 /2 (z - z-)(z - z+) 
with z± being the solutions of z\ — 2 (1 + B) z± + 1 = and 



(74) 



3 

B = £ (1 - cos(^)) + (aA) 2 /2 (75) 
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The non- vanishing terms as a goes to zero are Ci j2 (p = 0), which contain a linear divergence, 
and the first derivatives of Ci j2 {p) with respect to ps and p^ at p = 0, which are logarithmically 
divergent. In either case, the calculation reduces to the computation of the contour integral over z 
which can easily be performed taking into account that only the pole z = z~ lies in the unit circle. 
Furthermore, the e-prescription tells us that the pole of the quark propagator does not contribute to 
the contour integral when the Cauchy's theorem is used. The final result is, for the non-derivative 
contribution, 



d(2m + l,p = 0) 

Ci(2m,p = 0) = — (-) C F vl m 



by parity 
1 /a. 



7T 



Si (20) (m-1) + Si (11) (m-1) 



C 2 (2m + l,p = 0) = byparity 



C 2 (2m,p = 0) 



I C F „*"_L [Cs( 10 )(m-1) + Cs( 01 )(m-1) 



(76) 



The derivative with respect to ps at p = gives 

'**(*».»>•. _ „ byparity 



dps 

dd(2m -l,p) 
dps 
dC 2 (2m,p) 
dps 

0C 2 (2ro + l,p) 
dps 



— I c f v 3 



2m- 1 



2Cs^(m-l) - Id( n )(m-1) 



= by parity 



2m + 1 



(m — 1) + 

- (2m - 1) {Cs (10) (m - 1) + Cs( 01 )(m - 1) + Cs^(m - 1)}] (77) 



Finally, the derivatives with respect to p^ at p = yield 
'dCi(2m + l,pY 



dpA 
dd(2m,p) 

dpA 

dC 2 (2m + l,p) 
dpi 
dC 2 (2m,p) 

dpA 



by parity 



n c F vi m — 

TT J b 3 4 m 

Si( 2 °)(m-1) + Si( n )(m-1) + ^±lsi( 21 )(m 



by parity 

'a 



C F vl m 



2m 



-Cs' n '(m 



„ - 4- " - <78> 
where Si^\m), Cs^im) and Id (a/3) (m) are three-dimensional integrals which analytical expres- 
sions and numerical values for several m can be found in appendix A and Table A.l respectively. 
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We are now in a position to giving the expression of the heavy-quark self-energy on the lattice 
at any order in the velocity V3. In fact, it can be written as 



E(p, v) 



1 

a 



C F 



i=0 



+ 



7T 



+ — 



a 

+ 0(a s a) 



Cf 
4 

Cf 
4 



1 

00 
E 



i=0 

00 



Li=0 



SO i 



2i 
2i 



2 i 



e 



(2i) 



(«) 
40 



J 41 



log(oA) } 



-,(2i+l) 
J 30 



*,(2i+l) 
J 31 



log(oA) } 



(^3^3) 



(79) 



where the constants Sq30 40 can easuv De obtained from the results in eq.([7|) to (|78|), Their 
analytical expressions and numerical values are listed in appendix A and Table A. 2. 

To finish this section, we wish to discuss in detail the values of constants E^ and s!^ which 
determine both the one-loop wave function renormalization of the heavy-quark with velocity V3 and 
the anomalous dimension of the operator (v ■ D). 

We begin with y!i™\ For m = 0, it is the static heavy-quark wave-function renormalization, 
first computed in ref . P6|| . We reproduce their result E4} = 4 (see eq.( [A.6| ) and Table A.l). If 

(2m) 

m > 0, Tr A1 is determined by the pole part of the sum of the derivatives with respect to P4 of 
diagrams C.l and C.2. The terms that contain a logarithmic divergence are Si^ 21 ^(m) and Cs^ 11 ^ 

(2m) 

which appear in EL through the combination 



m 



,(2m) 2m + 1 



J 41 



x 



Si (21), 



m) I pole 



2m 



-Cs( n )(m) | pole 



(80) 



2 v ' n 2 

In appendix B, we calculate the logarithmically divergent part of si( 21 V) and Cs^^^m). Inserting 
eqs.( p.lO ) and (B.ll) into (p0|), we obtain that these integrals conspire order by order in the velocity 
to yield a vanishing coefficient E^ 2 "^ for m > 0. In other words, the coefficient of the logarithm 
of the wave-function renormalization of a heavy-quark moving with velocity 1)3 is independent of 
the velocity, as it should be. Moreover, the anomalous dimension turns out to be equal to that of 
the static theory, which in turn is the same as the one in the continuum. Therefore, we can say 
that the lagrangian (p2[) preserves both the infrared and ultraviolet behaviour of the non-expanded 
theory order by order in the velocity. 

Similarly, the pole part of the sum of the derivatives with respect to p% of diagrams C.l and 
C.2, determine Eg 2 [ m+1 \ In this case, there is only a term containing a logarithmic divergence, 
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Cs( n )(m), which appear in Eg^" 1-1 ' through the combination 

^ - ^ + * { 2C^(m) - #)( m ) } | P oie - ^ CB^frn - 1) | po ie (81) 



Again, substituting the expressions for the pole parts given in eqs.( B.ll| ) and ( B.12|) into (^), we 



observe that the coefficients of the logarithms conspire order by order in the velocity to produce 
a vanishing E^ for m > 0. The only logarithmically divergent term left is that for m = 0. 
Note also that Sgj is equal to E^J . This fact it very important because it implies that the 
renormalization constant of the operator (v ■ D) is finite (see next section), as it should be since 
this operator is conserved in the static theory. Again, consistency with the non-expanded theory 
is explicitly shown. 

Putting all these things together, we have 

w , 1 (<*s\ C F ^ (v 3 \ 2 * (2i) 



+ 



+ 



a£\ C F 

it J 4 

Os\ C F 

IT J 4 



|{0 2is «}- 41 ^ A) 

f;(^) 2i Eg +1 »}-41og(aA) 



(P3V3) 



Li=0 

+ 0(a s a) (82) 

which is one of our most important results. Notice that up to 0(d|), the heavy-quark self-energy 
( |82| ) coincides (within an error of less than 1%) with the one computed in the previous section 
using a different integration method (see eq.(^)). This fact makes us think that our numerical 
calculation is correct. 

5.3 Wave function and mass renormalization 

Having obtained the heavy-quark self-energy, we want to define and compute its wave function 
renormalization Zq (defined by Q R = Zq 1 ^ 2 Q b ) and mass renormalization SM. 

In order to get the renormalization constants, we study the heavy-quark propagator near on- 
shell including order a s corrections 

iH[P,V3) = (1 - S 4 ) (i P4U4 ) + (1 - S 3 ) ( P3 v 3 ) - S + 0(p 2 ) (83) 



where 
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The analytical expressions for S4 and £3 are readily computable from eq.(|8z 

If we impose on-shell renormalization conditions along with the normalization of the velocity 



1, it is easy to check that up to order a s [16] 



5M = -S (85) 
Zq = 1 + v\ £4 — v\ S 3 



1 + (v) x - 4 log(aA) + vl e£) + vl £ (sf { s - +2) - s£+1) 



i=0 

Z" 1 = 1 + v\ S 4 - v\ S 3 

2i 



(86) 



= ' + (?) ^ s © * K f) - } 

i=0 

where Z„ is the renormalization of the heavy quark velocity, first introduced in ref.(l(|], defined by 

v* = Z- 1 v£ (88) 

As we will see in the next section, Z v is a lattice effect that originates from the fact that the 
wave function and mass renormalizations are not sufficient to match the lattice and continuum 
amplitudes. In this sense, this 'velocity' renormalization can be interpreted also as the matching 
constant necessary to reproduce the physical amplitudes in the continuum from the ones on the 
lattice. We will return to this subject in the next section. 

For future use, it is convenient to write Zq and Z v as 

- -4 log(aA) + vl 

8=0 



z Q . 1 + ^)^[_ 41og(o „ + „Jgg) 2 - 2 f. 



(89) 



1/4 —\V t 

and the reader can find the numerical values of Zq 1 ^ in Table A. 2. 

To finish this section we wish to briefly comment on some particularly interesting characteristics 
of Zq and Z v . We observe that Zq is logarithmically divergent so that the corresponding anomalous 
dimension coincides with the one in the continuum and is independent of the velocity. On the other 
hand, Z v is finite because the logarithmically divergent terms coming from the self-energy cancel 
out exactly in eq.(^) order by order in the velocity. Finally, since lattice regularization breaks the 
(9(4) symmetry, the renormalization constants of the effective theory depend on the velocity of the 
heavy quark. This does not happen in a Lorentz invariant theory because there all quantities must 
depend on Lorentz invariants. 
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5.4 Renormalization and matching of (v ■ D) 

The 'velocity' operator (y ■ D) does not renormalize multiplicatively and, in general, we need to per- 
form subtractions of terms that diverge as powers of the ultra-violet cut-off a -1 . Specifically, we will 
demonstrate that (v ■ D) mixes under renormalization with the operator K(x) = Q\x)D^(x)Q(x) 
through a coefficient free of power divergences and with the operator l(x) = (x)Q(x) whose 
coefficient diverges as 1/a. 

It is convenient to proceed order by order in the velocity because in this way the matching can 
be understood better. 

Consider a single insertion of (v ■ D). The vanishing of the corresponding amplitude at zero 
external momentum (eqs,(|42|) and (|47l)) imply that (v ■ D) does not mix with the operator l(x) 
with a linearly divergent coefficient (i.e. proportional to 1/a). There is only a multiplicative 
renormalization of (v-D). In order to obtain it, we need the wave-function renormalization constant 
Z of the heavy quark in the static theory. It has been computed by many authors [Efl 



where = 24.48 (see eq.(82)). Then, the one- loop matrix element of the bare operator (v ■ D) 



between heavy quark states is given by 

\ v 3 \dp 3 J v 3 \dp 3 J \dpiJ o J 

<Xs\ C F (0) (1) 



= - [ 1 + { -f) -f ™ - S 30 ] J (t*Ps) + ■ • • (92) 

The dots indicate terms which vanish as p 2 for p — > 0, and therefore do not contribute to the 
on-shell renormalization. As anticipated in the previous section, the ultraviolet divergence of the 
vertex correction cancels the one of the field renormalization constant Z, leaving a finite term. This 
occurs because the operator (v • D) is conserved in the static theory. 

This result can easily be generalized for an odd number 2m+l of insertions of the operator (v-D). 
In fact, from (76) we learn that it does not mix with the operator 1 because the corresponding 



amplitude vanishes at p = due to the spatial parity invariance of the theory. Therefore, two-quark 
matrix elements of (v ■ D) do not contain linearly divergent terms proportional to 1. The same 
reasoning applies to the mixing with the operator K. Thus, (v ■ D) renormalizes multiplicatively 
with a finite renormalization constant at any odd order in the velocity. 
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Let us consider now the renormalization of the double insertion of (v ■ D). The amplitudes at 
zero external momentum are now non- vanishing, 



B 3 (0) + B 4 (0) + 0, 



(93) 



implying that there is a mixing of the double insertion of (v ■ D) with the operator l(x) with a 
linearly divergent coefficient. There is also a mixing of the double insertion of (v ■ D) with the 
operator K(x), because 



fdB 3 



+ 



8Ba 



(94) 



\dp 4 J ' \dp 4 / 

The mixing is finite because the logarithmic divergences of (dB^/dp^o and (dB^/dp^o cancel each 
other. This is true at any order in the velocity, as we demonstrated before. 

We have therefore the one-loop result 



Cp 

4 



Cp 

4 



(95) 



where the dots indicate terms which do not contribute to the on-shell renormalization. 



As before, these results can be extended to any order in the velocity by means of eqs.(|7q) and 
ffq) . We only give the final result for the one-loop renormalized lattice operator (v ■ D) 



[(v-D)] 



(i) 

Latt 



1 + 



1 / Ois 

a \ ir 



+ - — 



7T J 4 

C F 



Cf 2 



E 

i=0 



2 1 



+ 



7T 



4 



oo 

E 

i=0 



V- (2i) 



/ v (2i+l 
11(0) 



„(2i) 
J 40 



^3\ 2i y (2i) 



(96) 



where the superscripts (0) denote bare operators. 

To proceed further, we match the amplitudes of the bare lattice SHQET onto the amplitudes 



of the MS'-renormalized HQET (i.e. non expanded in v 3 ). To do this we need to know the two- 
quark amplitude of the operator (y ■ D) for an external momentum configuration on-shell. By 
direct computation, we observe that order by order in the velocity ^3 the sum of the relevant loop 
diagrams with insertions of ^-vertices vanish. Of course, there is a physical reason for this to 
happen: the operator (v ■ D) is conserved in the continuum static theory, therefore it does not get 
renormalized by interactions with gluons. Thus, we can write 



■d)Wg= [(v ■ D)] {0) + 0( 



(97) 
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In other words, the one-loop wave function renormalization of a heavy quark is independent of its 
velocity due to the fact that dimensional regularization is a covariant regularization. 

We can now perform the continuum-lattice matching by computing the ratio of the continuum 
amplitude to the lattice one. Prom eq.(p6|) and ( |97| ) we have that the physical operator (v ■ D) is 
related to the lattice bare one by 



s(t)> 



-,(2i+l) 
J 30 



,(2i) 
J 40 



1 /a,\ C 



a \ tt 

Os\ C F 
TT ) 4 



V4 



E 

i=0 



[l]Lott 



V- l^3\ V (M) 



[Kl 



Latt 



= z.„ 



[(v-D)] 



ci 



Latt 



\Latt 



c 2 [k; 



Latt 



[(v-D)] Latt 



(98) 



with obvious notation. As anticipated, we learn from the previous equation that Z v can be inter- 
preted as the lattice-continuum matching constant of the operator (v-D). 



An equivalent way of performing the matching is to expand the MS-renormalized propagator 
in the HQET and compare it order by order in the v 3 with the propagator in the SHEQT on the 
lattice. In the continuum we have 



iS(k) 



7 

"MS 



n MS 



iv^ki + v 3 k 3 



+ 



iki + e iki + e 



z ms , , ,y Z MS 
(-v 3 h) 

iki + e 



+ 



Z 



MS 



-v 3 k 3 ) 



1 



iki + e 



-v 3 k 3 ) 



Z 



MS 



+ 



Z MS —iv\ 



J MS 



ik& + e ik^ + e 2 



i&4 + e 



+ 



where Zj^g is the heavy quark field renormalization constant 



'MS 



I -l- (^) ^21ogWA) 2 



(99) 



(100) 



The bare lattice propagator in the SHQET is instead given (near the mass-shell) by: 
Z yZ , , ^ 1 _ /a.,\ Cf r „rni „m, \ yZ 



iS(k) 
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+ 



+ 



i&4 + e ik^ + e 



-v 3 k 3 



'oA Cf [s (o) _ I 



2/04 + e 



(-»***) 1 + ) 



Os^l CV rv (o) 
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ik^ + e 
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(101) 
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where Z is the field renormalization constant of the static lattice theory given in eq. (^Tj) (we omit 
for simplicity the mass renormalization). 

Matching at lowest order in V3 (static approximation) is realized by introducing a matching constant 
C of the bare lattice regulated field onto the MS renormalized field 

Qms = C Ql (102) 

where 

C = %2 = 1 + (^) % [21ogM 2 - £$] (103) 
At order v$, we must introduce a matching constant Z v defined by 

(v-D)j^ = Z v (v-D) BL (104) 
The comparison of eqs.(^) and ( |101| ) gives 

Z v = I+^^PS-Sg] (105) 

Matching at order requires to subtract from the double insertion of (v ■ D), the contribution 
proportional to K and the one proportional to 1, since they are absent in the HQET propagator. 
This means that there is a mixing of these operators in the lattice-continuum matching. Performing 
the subtraction above, we reproduce (|95|), 

This procedure can be iterated to higher orders in the velocity V3 leading to eq.(p7[). 

6 Renormalization of the heavy quark current 

In this section, we deal with the renormalization of the heavy-quark current 

J[x) = Q^(x,v)TQ(x,v') (106) 

describing the transition of a heavy quark with velocity v into a heavy quark of velocity v' . T 
stands for any of the 16 Dirac matrices. We specialize our computation to the most interesting case 
v' = (1,0) and v = (0,0,U3,^/l + u|). 

We will demonstrate that the weak current J(x) renormalizes multiplicatively with a coefficient 
that is only logarithmically divergent. We will show explicitly that the one-loop anomalous di- 
mension of J(x) depends on the velocity and coincides order by order in the velocity with the one 
computed within HQET in the continuum. 
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For the sake of clarity, we divide this section in two parts. In the first one, we compute the 
on-shell lattice matrix element of J(x) between heavy quark states up to order 0(v^). The second 
subsection is devoted to extend the previous result to all orders in the velocity. It is there where we 
will re-obtain the velocity-dependent one-loop anomalous dimension by summing all the diagrams 
with insertions of the operator (v ■ D). 

6.1 Matrix element of the current up to 

We start by considering the renormalization of the weak current J{x) with one insertion of the 
operator (v ■ D). We study the Green's function 

G(x, z) = J d 4 yd 4 w (0 | T [Q(x) Q(y)(v • D)Qt(y) J(w) Qt( z )] | 0) (107) 

The only non-vanishing diagram involved is drawn in Fig. 5. The amplitude of diagram D.l is given 
by 

9 2 C F v 3 Jq4 _ 2ip4 f d 4 k e-^ S m(h+P3) 1 



where p is the final momentum of Q and q is the (incoming) momentum of J. There is a potential 
logarithmic divergence, which is absent because the integral vanishes by parity at zero external 
momenta 

D x (p = 0,q = 0) = (109) 

Therefore, there is not any additional renormalization of the Green's functions of the form ( |107j ). 

Consider now the renormalization of the Green functions of J(x) containing a double insertion 
of {v ■ D). The relevant diagram is given in Fig. 5. The amplitude for D.2 is given, at zero external 
momenta, by 

_ g 2 C F vl f d 4 k r,(k) e~ ik * 1 
2 6 J (2vr) 4 (l-e-^ +e )4A 1 (A;) K > 

With the technique described in sec. 5.1, the previous integral is transformed into 
D2 = 1^"^ ( ~ 367T 1 / Ai{k) 3 ^ COS ^ 4 + 2 cos 2ki + cos 3fc4 ] 

rj(k) \ 



I d 4 k A 7?( /; ;j 3/2 + 3cosfc 4 + 2cos2A:4 ]- — / d 3 k 

J AilKI^ 127T J 



36^ 2 J Ax(fc) 2L ' " J 12tt7 A!(0,k) 2 J 

The logarithmic singularity of the amplitude is entirely contained in the first integral. The numerical 
computation gives 

a s \ Cf 2 r 2 i t \\i 



D2 = ( -f J -f vi [ - log(a\y + d 2 } (112) 
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where d 2 = -5.022. 

Adding to ( |112| ) the contribution from the external wave- function renormalization Zq , we obtain 
that the matrix element of the current J(x) between on-shell heavy-quark states is 

(c,v 3 | J | 6,0) = 1 - 9L [Z f + v\zf - \v\ log(aA)] 

+ 0(a s a,vl) (113) 

where Z^ = —19.95 is the old result for a static heavy quark and Z^ = —4.653. The reason for 
the introduction of the constants Z^ will be apparent in section 6.3. 

Note that the logarithmic divergence from the vertex diagram for two static heavy quarks (lower 
order in V3) exactly cancels the one from the external wave- function renormalization resulting in 
a finite lowest order correction to the matrix element of the current zj?' . The physical reason 
for this to happen is that the flavour conserving current, i.e. the current J(x) for equal velocities 
v = v' or equivalently v% = 0, is conserved in the HQET and so its anomalous dimension must be 
zero. Therefore, the anomalous dimension of the current J(x) starts from in an expansion in 
the velocity. 



6.2 Matrix element of the current beyond 0(vf) 

The only non- vanishing diagram we need to calculate now has the same structure as those in Fig. 5 
but with m insertions of the operator (v-D). We will denote it by E. All other possible one-particle 
irreducible diagrams vanish due to parity. Again, we do not consider the insertion of ^-vertices 
because the net effect of all such vertices is to multiply the original diagram by l/vj 1 , with m the 
number of i^-vertices. The demonstration of this assertion is similar to the case of the self-energy 
and so we do not repeat it here. 

The amplitude corresponding to the diagram E is 

d 3 k 



dz z —z 



(114) 



2vri z [1 + e - e~ ipi z] m+2 (z-Z-)(z- z+) 
where z = e iki and the contour integral is along the unit circle. 

The non- vanishing term as a goes to zero is E(p = 0), which contains a linear divergence. 
First derivatives with respect to the external momentum give rise to terms of order 0{a) that do 
not contribute to the on-shell renormalization. The computation reduces to the evaluation of the 
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contour integral over z which can be easily performed taking into account that only the pole z = Z- 
lies in the unit circle. Furthermore, the e-prescription tells us that the pole of the quark propagator 
does not contribute to the contour integral when the Cauchy's theorem is used. The result is, 

E(2m + l,p = 0) = by parity 



E(2m,p = 0) = - — C F v 



,,2m 



it J ^ r " 3 4"i+l 
Si (20) (m - 1) + Si (11) (m - 1) + Si (21) (m - 1)1 (115) 

For m = 0, the static case, the amplitude greatly simplifies 

E(0,p = 0) = -^^ld^(0) (116) 

The integrals Si (a/3) (m) and ld (11) (0) are defined and their numerical values tabulated in appendix 
A. 

The matrix element of the current J(x) between heavy-quark states is logarithmically divergent 
for m = and also for m > because so is Si^ 20 ^ (to— 1). The coefficients of the logarithms can easily 
be extracted order by order in the velocity from eqs.( [B.10| ) and ( |B.12j ). Adding the contribution 
from the wave-function renormalization of the external states (see eq.(|86l) and (^)), we have that 
the matrix element of the current can be written as 



< c ,„ 3 uM> = i + (^)f g{-tf*f° 

— ) 2 



V3\ 
V 4 



J 



1 r 

Oifl 



(2i + 1) 



log(aA) 2 \ + 0(a s a) (117) 



where can easily be evaluated from eqs.( |115[ ), (116) and (|89|). In order to simplify the com- 



putation in section 6.4 of the relation between the Isgur-Wise function on the lattice and in the 



continuum MS, we have expanded the finite contributions in powers of v% instead of v^/v^. This 
is achieved by noting that 

-) 2m .f; ( _y(i±!^ f «> (118) 

(n) 

For the numerical values of constants ' we refer the reader to Table A. 2. 

The only subtlety in this calculation is the fact that the wave function renormalization at lowest 
order in the velocity (i.e. the static case) contributes to the matrix element with a coefficient twice 
the one of higher velocity orders, which is 1/2. The reason is that we consider the b quark static 
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and the c quark moving with a small velocity v 3 . Therefore, only the latter, as a consequence of 
its interaction with the gluon field, gets a velocity dependent wave function renormalization which 
lowest order is the corresponding to a static heavy quark. 

The interesting thing is that the infinite sum in front of log(aA) can be evaluated simply recalling 
that 



2u 



1 



u 



,2i 



( 2 * + 



1) 



Substituting (|TTg) into (glT]), we get 

(c, v 3 I J I 6, ) 

+ 0(a s a 



1 ; 



-2 + V -± log 

^3 



^4 + ^3 

v 4 - v 3 



log(aA) 2 - J2 { 



v 2i Z ( 2i ) 
V 3 Zj^ 



i=Q 



(119) 
} 

(120) 



It should be stressed that the infrared structure of the matrix element of the heavy-heavy current is 



the same as the one evaluated within the continuum Georgi's theory in ref.[27], as it should be. In 
other words, the expansion in the velocity reproduces the correct infrared behaviour of the theory 
once we sum all orders in the velocity. This is a check of the consistency of our approach. 

(2) 

Another check is provided by the comparison of the numerical value of the constant listed 
in Table A. 2 (computed by direct contour integration) and the one given below eq.( |113| ) (computed 
by integration by parts). They coincide within error bars. 



6.3 Lattice-continuum matching 



In order to match the two-quark amplitude of the current to its counterpart in the continuum, we 



need to compute the matrix element of the current in the HQET renormalized in the MS scheme. 



This has already been done by the authors of ref.[27] and here we only quote their final result 



(c,v 3 \J\b,0) Jgs =l + _± _jl 2-^lo. 



a. 



7T J 4 

+ 0(a 2 s 



V3 



V4 + V3 



:i21) 



where \x is the renormalization point. 



By forming the ratio of (121) to ([L20|) we get the factor that relates the matrix elements of the 



heavy-quark current J(x) in the lattice and in the continuum MS renormalization schemes 



(c, v 3 I J I b, 

(c, v 3 I J I b, ) LaU 



1 ! - < v) X 



-2 H log 

V 3 \V4-V 3 



log(a^) 5 



+ 



(JXJ 

£{^f } } 



+ 0(a 2 s 



(122) 



i=0 
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As expected, the infrared regulator A disappears in the matching because the lattice theory and 
its counterpart in the continuum have the same infrared behaviour. It is the ultraviolet one that is 
different and ( |122[) takes this discrepancy into account. 

Let us discuss now the on-shell renormalization of the lattice SHQET in the real space ]28|] 
instead of momentum space p6|| as we have done up to now. These renormalization schemes differ 



on the lattice and the relation between them has been clarified in [29|. A clear exposition can be 
found in which we will follow almost verbatim here. 

Consider the bare propagator of the heavy quark moving on the lattice with a velocity v 3 along 
the z-axis as a function of time and momentum p. We will call it iH(t,p). This propagator can 
be obtained by performing the Fourier transform with respect to P4, the fourth component of the 



external momentum, of the propagator in the momentum space (82). For large euclidean time and 
in the continuum limit, iH(t,p) reduces to 

iH(t,p) = Zq 6 -^- exp[-(t + l)/vf (SM + v§ ■ p 3 ) } (123) 

Note that for the momentum space renormalization conditions ( |85|) to (pf\), Zq is multiplied by an 
exponential with (t + 1) instead of t in the heavy quark propagator. 

On the other hand, in lattice simulations one fits the correlation functions to an exponential 
evolution in euclidean time with t instead of (i+1) (real space renormalization scheme). Therefore, if 
we do not modify the momentum space renormalization conditions appropriately, we will not take 
into account the correct wave function renormalization giving rise to a wrong lattice-continuum 
matching. The solution is to take a shifted wave function renormalization Zq related to the old 



one by |26| 

Z Q = Zq - — (124) 

which tells us that he discrepancy between the momentum and the real space renormalization 
schemes is finite and is given by the mass renormalization. In Table A. 2, we have tabulated the 
values of the renormalization constants Zq and Z^ in the real space lattice scheme. 



7 The Isgur-Wise function 

In this section we determine the relation between the value of the derivatives of the Isgur-Wise 
function at the zero recoil point, ^ n ^(l), measured on the lattice and its physical counterparts in 
the continuum MS. 



28 



Isgur-Wise 

LLtJI 1 Vex LI V Co 


Numerical Coefficients 


1 VLatW SLatUJ S^atv 1 / 


Ae (1) (« _1 ) 


-9.31 -19.95 0.0 0.0 
3.44 -18.62 -19.95 0.0 
88 20 10 31 —27 92 —19 95 


Af\a^) 

AfV 1 ) 
Al^ia' 1 ) 


-24.88 0.0 0.0 0.0 
17.35 -49.76 0.0 0.0 
10.56 52.04 -74.64 0.0 



Table 1: Numerical values of the constants determining the continuum-lattice matching of the first 
derivatives of the Isgur-Wise function. A factor multiplying all entries is understood. 



In fact, taking fj, = a in ( 122 ), we find 



i + 



a s (a i )\ C F 



£,Lat{Vi) \ vr 



(125) 



we have 



where the matching constants are defined in eq.( [ll7[) . Note that by setting fi = 
got ride of the logarithms that appear in the lattice-continuum matching. 

Substituting the expansions of both the Isgur-Wise function in the continuum and the one on 



the lattice in powers of V3 into ( 125 ) and demanding consistency order by order in V3, we get 

/ a s (a _1 )\ Cf 



A^ia- 1 ) 
A^ 2 \a- 1 ) 

A^\a~ l ) 



7T 



{^ o) d\(i) + 2^ 2) } 



a s (a 



/ ^{^ ) a(l) + 4Zf^) ( l) 

\zf) 

)^{^ ) ^l) + (6Zf + 24^)eE(l) 



7T 



+ 2zf + 



a, (a 



7T 



+ 62f + 24^+48 Z$>} 



(126) 



(127) 



(128) 



with £( n )(l) being the nth derivative of the Isgur-Wise function with respect to V4 at the zero recoil 



point v A = 1 and A^)( Mo ) = l„ ^Zi 1 )- 



; (n) 



MS' 



Eqs.( |l26"|) to (|128[) are our most important results. They give the one-loop relation at the 
scale /x = o -1 between the lattice measures of the derivatives of the Isgur-Wise function and their 
physical values. It should be stressed that equivalent expressions can be obtain for the real space 
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Isgur-Wise 
derivatives 


R. G. 


correction 




a' 1 = 2 GeV a" 1 


= 4 GeV a' 1 


= 6 GeV 


Azjp(m) 


0.668 


2.111 


3.028 


Azf\m) 


-0.248 


-0.789 


-1.135 


Azf\m) 


0.155 


0.489 


0.701 



Table 2: Renormalization group (R. G. ) corrections to the constants determining the contin- 
uum-lattice matching of the first derivatives of the Isgur-Wise function at the scale fi = fn for 
several lattice spacings a. The two-loop anomalous dimension of the heavy-heavy quark current 
has been properly included. A factor (^) —g- multiplying all entries is understood. 



renormalization scheme by replacing zj^ by Z^\ In Table 1, we give the numerical values of the 
coefficients of £^(1) in A^ n \a~ l ) both for the momentum and real space (denoted with a bar) 
renormalization schemes. 

As we mention in Sec. 4, the values of the matching coefficients in Table 1 depend on the contin- 
uum renormalization scheme. Although we expect their numerical values not to change very much 
in a different renormalization scheme, consistency requires to properly include the contribution of 
the two-loop anomalous dimension. In addition, we give the renormalization group evolution of the 
derivatives of the Isgur-Wise function from the scale a -1 to a generic renormalization point //. 

The one-loop anomalous dimension 71 of the heavy-heavy quark current is velocity dependent 
and has been calculated in eq. (|l20| ). On the other hand, the two- loop anomalous dimension 72 has 
been computed in ref . 1 30 1 in the MS scheme. Expanding both 71 and 72 as a power series in v% 
and inserting the result in eq.(|35"|), we find that the renormalization group corrections AZ^ n \p) to 
the matching constants zj^ are 

/ a s (a _1 ) 



7T 



a s ( a 



7T 



^AZf(M) 



+ 



a s (a~ 1 ) a s (n) 



7T 



7T 



2vr 2 + 



18 



N F 



29 02 
6 Px 



(129) 



~^~lh { ~^ lo s(«s(>)/«s( a X )) 



7T 



7T 



4vr2 
5 



N F 



61 

40 



5 (3t 



(130) 



30 



/^.(a" 1 ) a s (^)\ [tt 2 , 4 Ar 43501 , 8 ft"L. }1 . 



7T 7T 



vr 2 4 AT 43501 8 02 

I Nf 1 — 

35 63 37800 35 ft 



In Table 2, we report the numerical renormalization group corrections to the constants zj^ for 
several values of a -1 at the physical meaningful scale [i = rn where fn = mBTnr,lijriB ~ 1-4 

GeV, i.e. the reduced mass of the B- and D-meson. The number of active quarks is 3 because both 
the b and the c quarks are taken to be static sources of color. The values in Table 2 must be added 
to those of Table A.l to obtain the matching constants at the scale fn. 



8 Power divergences 

In this section we prove that the renormalization of £^(1) is not affected by ultraviolet power 
divergences. The argument presented does not rely on any perturbative expansion and is based 
only on the symmetries of the lattice SHQET. 

A given operator O can mix with lower and equal dimensional operators O' allowed by the 
symmetries of the (regulated) theory. The mixing coefficients are proportional to the appropriate 
power of the ultraviolet cut-off l/a to account for the dimension. If the dimension of O and of 
O' are the same, the mixing coefficients contain in general logarithmic divergences, of the form 
log a. The computation of £^(1) involves single and double insertions of the velocity operator 
(v ■ D) = Q\v ■ D)Q. The only possible linear divergence in (v ■ D) is through the mixing with the 
operator 1 = Q^Q 

k 

< (v ■ D) > = — < 1 > + at most logarithmic terms (132) 

where k is a coefficient which has a perturbative expansion in a s . Such a mixing is however 
impossible because of the spatial parity invariance of the theory. Then, we have k = 0. 

The double insertion of (v ■ D) can also mix with 1, and in this case the mixing is not forbidden 
by any symmetry 

_> _> c 

< (v ■ D) (v ■ D) > = - < 1 > + at most logarithmic terms (133) 

where now c / in general. In sec. 5 we have checked this result with an explicit one-loop compu- 
tation. 
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Therefore, the two- and three-point correlators defined in eqs.(|l4|) and ( |T6| ) respectively, can be 
written as 



C^C*' - *) = - {t' -t) C$\t' -t) + at most logarithmic terms (134) 



C 3 (t,t') = — {t' — t) C\ (t,t') + at most logarithmic terms 
c 



a 



where the superscript indicates the order in the velocity. The derivative of the Isgur-Wise function 
is given by the following combination of correlation functions (see eq.(|l9|)) 

Substituting eqs.( |134[ ), we see that the linear divergence cancel in the expression for the £ (1) (1), 



as anticipated. We notice that the argument given is non-perturbative, and it is confirmed by the 
explicit one-loop computations of the previous sections. 

This argument, however, does not hold for higher derivatives of the Isgur-Wise function. In 
fact, consider for example, the second derivative of this function with respect to v^. In this case, 
we must deal with four insertions of the operator (v • D). Therefore, the correlation functions C^j 
will contain a linearly divergent contribution just as C^. From eq.(^) we see that the correlation 
functions C^l enter the expression for £^(1) through the combination 

C W a® - Of C« (136) 

which again is at most logarithmically divergent because the poles 1 /a cancel out. There is, however, 
a second contribution to £^(1) that is not free from linear divergences, namely, 

(C 2 (2) ) 2 cf -Cfcf Cf (137) 



Substituting (134) into (137) we obtain that the term 1/a 2 cancels out in this combination but 
that proportional to 1/a survives giving rise to a linear divergence as a goes to zero. Therefore, 
the computation of £( 2 )(1) requires the subtraction both from and of a linear divergence 
as in the case of the self-energy of a quark. 

The same reasoning can be applied to higher derivatives of the Isgur-Wise function. As we 
increase the order of the derivative, the power of the divergence also increases and thus non- 
perturbative subtractions from the correlation functions of lower velocity degree are necessary to 
obtain reliable results from a numerical computation on the lattice. 
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9 Conclusions 



We have studied the lattice renormalization of the effective theory for slow heavy quarks, which 
allows to compute the slope of the Isgur-Wise function at the normalization point, with 
Montecarlo simulations. We showed that the lattice-continuum renormalization constant of £^(1) 
does not contain any linear ultraviolet divergence, but only a logarithmic one. This implies that 
the matching of £W(1) can be done in perturbation theory and it is not necessary to perform any 
non-perturbative subtraction. The lattice computation of the slope of the Isgur-Wise function using 
the effective theory for slow heavy quarks is therefore feasible in principle. 

The one-loop lattice renormalization constants of the slow heavy quark effective theory have 
been computed to order v 2 together with the matching constant of ^^(l), which relates the value 
of this form factor measured on the lattice to its physical counterpart in the continuum. 

We have demonstrated that the effective theory for slow heavy quarks reproduces the infrared 
behaviour of the original (non-expanded) theory order by order in the velocity. This means that 
we are dealing with a consistent expansion of the effective theory for heavy quarks. 

We have analysed the lattice effective theory for slow heavy quarks also to higher orders in the 
velocity. Unfortunately, the renormalization of the higher derivatives of the Isgur-Wise function, 
£( n )(l) for n > 1, is affected by ultraviolet power divergences. The lattice-continuum matching 
of £( n )(l) is therefore much more involved than in the case of ^^(l). We stress however that 
the higher derivatives of the Isgur-Wise function are much less important than the first one. The 
renormalization problems of the slow heavy quark effective theory which arise in higher orders do 
not constitute therefore a serious limitation for its phenomenological applications. 

We hope that the results of our analysis may encourage the scientific community to carry out 
the numerical simulation of £^(1) using the slow heavy quark effective theory. We believe that 
this theory can be a source of interesting physical results. 
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Table captions 



Table 1: Numerical values of the constants determining the continuum-lattice matching of the 
first derivatives of the Isgur-Wise function. 

Table 2: Renormalization group (R. G. ) corrections to the constants determining the continuum- 
lattice matching of the first derivatives of the Isgur-Wise function at the scale \i = m for several 
lattice spacings a. The two-loop anomalous dimension of the heavy-heavy quark current has been 
properly included. A factor (jjf^j multiplying all entries is understood. 

Table A.l: Numerical values of three-dimensional integrals for several values of m, the order in 
the expansion in powers of the velocity i>3. 

Table A. 2: Numerical values of the constants entering the continuum-lattice matching of the 
heavy-quark current for several values of m, the order in the velocity expansion. 



Figure captions 

Figure 1: The diagrams contributing to the one-loop heavy quark self-energy with one insertion 
of (v ■ D), denoted by a crossed circle. The double line represents the heavy quark with velocity 
(0,0, t*). 

Figure 2: The diagrams contributing to the one-loop heavy quark self-energy with two insertions 
of (v ■ D), denoted by a crossed circle. The double line represents the heavy quark with velocity 
(0,0, t*). 

Figure 3: The non-vanishing diagrams contributing to the one-loop heavy quark self-energy 
with m insertions of (v-D), denoted by a crossed circle. The double line represents the heavy quark 
with velocity (0, 0, vs). 

Figure 4: The non-vanishing diagrams contributing to the one-loop heavy quark self-energy 
with m insertions of (v-D), denoted by a crossed circle, and n insertions of the operator D^, 
represented by a full circle. The double line represents the heavy quark with velocity (0,0,^3). 
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Figure 5: The non-vanishing one-particle irreducible diagrams contributing to the one-loop 
vertex of the current J(x) with one and two insertions of (v ■ D), denoted by a crossed circle. The 
double line represents the heavy quark with velocity (0, 0, V3). The full line stands for a static heavy 
quark. 
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Appendix A Analytical expressions and numerical values of 

loop integrals 



The renormalization of both the heavy-quark current and the SHQET lagrangian can be written 
in terms of a few three-dimensional one-loop integrals. In this appendix, we give their analytical 
expressions and numerical values up to O(uf). 

We define 

,2/ 



SiM)( m ) = — r d 3 k sin {h) ^my m 

2vr J- n B a A( 2 + B)B} 3 



Cs^)(m) = 1 f ™Ah/2) g(p)m 



2vr 



B a \ [{2 + B)Bf 



Id (a/3) (m) = /* A ; 1 = 5(ff) m (A.l) 



B a J[(2 + B)B]P 



where for the present calculation a = 0, 2 and /3 = 0, 1, and 



sin 2 (fc 3 )(2 + ff) 
=(jB) = B 



B 

1 + 



2 



4 sin 2 (A;. 



3 J 



^2 



(A.2) 



^{2 + B)B 

with z_ the solution of z 2 _ — 2 (1 + B) z~ + 1 = with | z_ |< 1. Note that the function H(2?) 
is infrared convergent. This fact will be used in appendix B to subtract the infrared divergent 
behaviour from the integrals in eq.( |A.l| ). 

Other integrals as, for example, the one with a factor cos(fes) instead of cos 2 (&3/2) can trivially 
be reduced to linear combinations of the integrals defined in eq.( |A.l| ). 

Obviously, these integrals must be evaluated numerically. However, care should be taken when 
infrared divergences appear as in Si*- 21 ^(m), Cs^^m) and Id^^m). In fact, in this case we 
cannot take A = 0. The logarithmic infrared divergence must be subtracted before computing 
them numerically. This has been done for arbitrary m in appendix B. We refer the reader to this 
appendix for details. Other integrals are infrared finite for any value of m and thus can safely be 
computed by means of, for example, a Monte Carlo routine. 

Now, we give the analytical expressions of the heavy-quark self-energy (see eq.(|82[)) 

1 



V>(2m) 



• ' a m— 1 



Cs ( 10 )(m-1) + Cs( 01 ) (m — 1) 

(A.3) 



isi( 2 °)(m-l) - lsi( n )(m-l) 
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Velocity Power 


Infr. Diver. 


Type 


B-Factor 


Tfl — 


Tfl — 1 


m — 2 


log(aA) 


Ol 




8.284 


26.148 


178.85 


n 
u 


Si 


(11) 


3.367 


15.20 


78.8 





Si 


(01) 


5.791 


24.877 


120.76 





Si 
Ol 


{ZU) 

\ ZL ) 


6.771 
— U.Uoo 


29.2 
y.oo 


148 

78 O 

( o.y 


n 

u 

2 ±m+l 

2m + 3 


Ci 
Ci 
Ci 


(10) 

(11) 
(01) 


2.485 
7.298 


041.0 

12.93 
21.21 


1 zLQ 

72.9 
88.0 


o 

4»Tt+l 
1 — T 

2m + 1 



Id 


(10) 


19.95 


39.2 


241 





Id 


(11) 


4.526 


20.21 


108.5 


2m + 1 


Id 


(01) 


12.23 


34.9 


152 






Table A.l: Numerical values of three-dimensional integrals for several values of m, the order in the 
expansion in powers of the velocity V3. 





Velocity Power 


Constant 


m = 


m = 1 


m = 2 


m = 3 


m = 4 


m = 5 


m = 6 


v (m) 


-19.95 


0.0 


15.57 


0.0 


8.20 


0.0 


7.75 


v (m) 


24.48 


0.0 


7.60 


0.0 


7.55 


0.0 


8.65 


v (m) 
^30 


0.0 


12.67 


0.0 


7.28 


0.0 


8.37 


0.0 


7W 


24.48 


0.0 


-5.07 


0.0 


0.27 


0.0 


0.28 


Z Q 


4.53 


0.0 


30.45 


0.0 


-7.10 


0.0 


-0.17 


7H 


11.81 


0.0 


0.32 


0.0 


-0.82 


0.0 






19.95 


0.0 


4.65 


0.0 


-1.59 


0.0 


-1.04 




0.0 


0.0 


12.44 


0.0 


-5.28 


0.0 


2.42 



Table A. 2: Numerical values of the constants entering the continuum-lattice matching of the 
heavy-quark current for several values of m, the order in the velocity expansion. 
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s (2m+l) 



1 



30 Am— 1 



2m + 1 



{Si (10) (m - 1) + Si (01) (m - 1) + Si (11) (m - 1)} 0(m) 



s (2m) 



(2m - 1) {Cs (10) (m - 1) + Cs (01) (m - 1) + Cs (11) (m - 1)} 6(m) 

2m ±l { 2 Cs( 11 )(m) -Id( n )(m)} + \ 
1 



+ fH^± {2Cs( n )(m) - Id( n )(m)} + ± Id( 01 ) (m) 5 mfi 
Si (20) (m- 1) + Si (11) (m- 1) 



(A.4) 



40 Am- 1 



+ 2m + 1 Si (21) (m - 1) - (2m - 1) Cs (11) (m - 1)] (A.5) 
which are supplemented with the old results for a static heavy quark 

4°) = -I d ( 10 )(0) = Id( 1Q )(0) + ld( n )(0) (A.6) 

In table A.l, we list the numerical values of the three-dimensional lattice regularized integrals 



( A.l ). These quantities has been evaluated using both a Monte Carlo and a lattice integration 
routine. Errors are at most O(l) in the last decimal place. 

In table A. 2, we present the numerical values for the heavy-quark self-energy, the wave function 
renormalization, the velocity renormalization and the lattice-continuum matching constants for the 
Isgur-Wise function. As before, errors are at most O(l) in the last decimal place. The constants 
with a bar are in the real space lattice renormalization scheme while the others has been calculated 
in the momentum space scheme. 

Appendix B Infrared subtraction of loop integrals 

Feynman integrals appearing in lattice perturbation theory must be evaluated numerically be- 
cause they are too much complicate to be integrated analytically. The trouble arises when these 
integrals are divergent as A, the fictitious gluon mass, goes to zero. In order to compute divergent 
lattice integrals, one has to subtract from the integrand an expression which has its same infrared 



behavior 23]. Doing so, the integral to be calculated can be expressed as a sum of an infrared 
finite one, in which we can safely set A = 0, and a second integral which contains the divergences 
of the original one. The former can be evaluated numerically while the later must be computed 
analytically to explicitly display the terms which depend on the infrared and ultraviolet regulators. 

In this calculation, all four-dimensional integrals can be reduced to three-dimensional ones by 
integrating over the zeroth component of the loop momentum or performing an integration by 
parts. This simplification is possible because of the simple structure of the heavy quark propagator 
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that only depends on the zeroth component of the momentum. The integrands of the resulting 
three-dimensional integrals turn out to be algebraic functions of B, defined by 

3 AV 



a=l 

Therefore, it is convenient to know the infrared limit of B itself and some other functions of it. For 
| k a |<S 1 we have 

B « i^ + A^l-l^ + A 2 )] 



[1 - -^(Jfe' + A 



2 U-^] 



V^T+S) 3 - 1 (£ 2 + A 2 ) 3 / 2 4 

WW+W - 1 (fc 2 + A 2 ) 5 /2 [1 " 12^" 3 A ' ] (R2) 
The previous expansions are almost all we need to extract the logarithmic infrared divergence of 
our one-loop integrals at every order in the velocity. 

As we saw in appendix A, the only infrared divergent integrals are Si^ 21 ^(m), Cs( n )(m) and 
Id( n )(m), defined in eq.( |A.l| ). In order to numerically compute these integrals, we find their 
infrared behaviour using ( p.2| ) and then construct the corresponding regularizing integrals which 
are 

si<-V)„ = < -t, /?'(F5^ ,( "'- , ' ) 

/it {.2m 
d 3 k ^ Qiir 2 _ k 2 ) 

/it 1.2 m 

,/ 3fc (P+ ^3/ 2+W ^ 2 - fc2 ) ( B - 3 ) 

where we perform the integration on a 3-sphere of radius ir to take advantage of spherical symmetry. 
Of course any other radius R > would be equally good. 

The three-dimensional integrals in ( |B.3j ) are much simpler to be evaluated than the original 
ones. In fact, the best thing we can do is to separate the radial and angular integrations expanding 
fa as a Gegenbauer series 

fc n _ fe n I>) n [ y? (n-2j + v) ~ 
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where v is related to the space dimension D by D = 2{y + 1), therefore v = 1/2 in our case. 
Inserting (|B,4j ) into ( [B.3D , we have 



Si (21) (m) //? 



(2m + 3) 
1 

(2m + 1) 
1 

(2m + 1) 



I R (m + 2,m + 2) 
lR(m + l,m + 1) 



Iii(m + l,m + 1) 



where I R (a,f3) is the following radial integral 



dk- 



k 



2n 



+ A 2 ) 3/2+/3-l 

The integration limits are a consequence of the Heviside function introduced in (|B. 

The radial integral I/j(m,m) can be evaluated by noting that 

A 2 d 



(B.5) 



(B.6) 



lR(m,m) = Ir^tu — l,m — 1) + 



along with 



(3/2 + m-2) d\ 2 
I R (1, 1) = log(27r) - 1 - log(aA) 



I R (m - l,m- 1) 



The result is 



m— 1 



I R (m,m) = log(2vr) - ^ 



1 



log(aA) 



n=0 ( 2 " + ^ 

Putting all the formulas together, we arrive at the following infrared subtracted basic integrals 



(B.7) 
(B.8) 
(B.9) 
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TO 
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2^ 
24 
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2(m+l) 
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Cs( n )| 



TO 



2vr 



cos 2 (fe 3 /2) 



B^{2 + B)B 



(2n + 1) 
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log(aA) 



(B.10) 



2fe 
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(^ 2 - A; 2 ) 



4771+1 
(2m + 1) 



log(2vr) - X] 
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(B.ll) 



+ 



4771+1 
(2m + 1) 



log(2vr) - ^ 



1 



(jfe2 + A 2 ) 3/2+m 

log(Aa) 



(B.12) 



n= (2- + 1) 

Using the previous equations, we have calculated the numerical values of these divergent integrals 
which are tabulated in Table A. 2. Moreover, the coefficients of the logarithmic divergence determine 
the anomalous dimension of the heavy-heavy quark current and the running of the derivatives of 
the Isgur-Wise function. 
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